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Abstract 

We review a few useful concepts about polarization measurements in the quantum domain. 
Using a perfectly general formalism, we show how to build the quantum counterpart of some 
classical quantities like Stokes parameters and Mueller matrices, which are well known in classical 
polarization-measurement theory. 
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The effective polarization density matrix 



In these notes, our first aim is to re-derive the formulas for the effective density matrix 
p c g we introduced in without using the concept of Stokes parameters and without using 
a circular polarization basis. To begin with, let us establish our notation. In order to make 
the results comparable with the ones in Refs. Q], we adopt a relativistic notation. All 
formulas are given in natural units (c = H — 1) and all quantized fields are transverse. In 
this context, single-photon plane-wave states are denoted by |k, À) where k is the spatial 
part of the four momentum k = (k°, k), k° = |k| = u, and À = 1, 2 is the linear polarization. 
They are created from the vacuum state |0) by the corresponding creation operator aj^(k) 
which, together with the annihilation operator à\(k), satisfies the commutation relation 

[a A (k),4,(k')] = (27r) 3 2fcW (3) (k-k'). (1) 

With this convention, the normalization condition between two plane-wave states reads 

(k,A|k',A'> = (2n) 3 2k°5 xy 5^(k-k'), (2) 

which can be obtained directly by calculating the expectation value with respect to the 
vacuum state of the left side of Eq. ((TJ). Then the resolution of the identity can be written 

as 

1 = |0)(0|+ / dk^2 |k, A> <k, A| + ^ {multiparticle states} , (3) 

J A=l 

where (ik denotes the Lorentz-invariant measure 

dk = - — — — — . (4) 

(2n) 3 2k° v ; 

The plane-wave states |k, A) are quite special since they are eigenstates of the linear mo- 
mentum of the field. More generally, a single-photon state can be described by the density 
operator 

1,2 f 

P = J2 ^kJk' P Av(k,k')|k,A)(k',A'|, (5) 

A,A' J 

where p\y (k, k') is a Hermitian positive semidefmite 2x2 matrix. While the linear mo- 
mentum and the polarization degrees of freedom are decoupled for photons in the plane- 
wave states |k, A) this is not the case for photons described by p which spans an oo 2 - 
dimensional vector space (2 polarization directions for each wave vector k). Therefore, if 



one want to deal with polarization degrees of freedom only, it becomes necessary to reduce 
p to an effective finite-dimensional representation p c g (2 x 2 or 3 x 3) Q]. Now we present a 
derivation of p e g which enlightens the geometrical aspects of the problem. To this end, we 
have to do some manipulations on the fields. In the Coulomb gauge, which is assumed to 
hold throughout this notes, the free field vector potential operator À(r, t) is written 

r 2 

À(r, t)= dk [ e(A) (k)a A (k)e- í(fe0í - k - r) + h.c] , (6) 

^ A=l 

where the real unit vectors e^(k) {linear polarization) satisfy the following orthogonality 
conditions: 

e( 1 )(k). e ( 2 )(k) = 0, e«(k)x e ( 2 )(k) = ^. (7) 

Since we are working with the transverse field (A°(r,t) = 0), we can assume an Euclidean 
mètric in the spatial 3-dimensional space and make no distinction between high and low 
indices. It is useful to define e( 3 )(k) = k/A; and to build the k-dependent complete basis 
£(k) = {eW(k), e( 2 )(k), e^ 3 ^(k)}. Let us consider now three unit vectors x, y, z of a Cartesian 
coordinate system or, more generally, a set 1Z of three real orthogonal unit vectors TZ = 

{ e (D )e (2) ?e (3)} ; 

e M. e ®=6 ab , {a,b= 1,...,3), (8) 
which form a complete basis in the ordinary Euclidean 3-dimensional space: 

£eM:eM = l * j^e^ef^S^, (ij = 1, . . . , 3), (9) 

a=l a=l 

where the unit dyadic 1 has been written as a sum of dyadic products : (a = 
1, 2, 3). The vectors define three orthogonal spatial orientations and they are independent 
from the momentum direction k/k . However, for a given k, one can write the orthogonal 
transformation A(k) between the two basis £(k) and 1Z as 

A„ t (k) = e("' ■ e">(k) 

= X>, w 4"(k). (10) 

i=l 



Then we can write 



e (A) 



(k) = ^ e WA 6A (k), (A= 1,2,3), (11) 



6=1 
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and insert this formula in Eq. (JüJ) in order to obtain: 

3 



À(r,í) = 5>Wi 6 (r,í), (12) 

6=1 

where we have defined the ò-th component of the field A(r, t) as 

r 2 

À b (r,t) = /Jk^[A 6A (k)a A (k) e -^° í - k · r )+h.c.] 

J A=l 

(13) 

where we have defined the transformed annihilation operators A b (k) (ò = 1, 2, 3), as: 
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A(k) = AbA(k)a A (k), (b = 1, 2, 3). (14) 

A=l 

It is easy to check that these operators satisfy the following commutation relations 

[À a (k),Àl(k')] = (2n) 3 2k°A aa ,ó (3 \k - k'), (15) 

where we have defined the transverse Kronecker symbol A aa / as 

k k i 

A aa / = 5 aa > - (16) 
As expected, the longitudinal part — ? of A aa /, spoils the canonica/commutation relation. 



In the quantum theory of photo-detection it is a Standard practice |3| to define the positive 
frequency operators V&(r, t) (b = 1,2,3) as 

V b (r,t) = J Jkv^A(k)e- i(fc0í - k - r) , (17) 

and such that Yll=i ^ft~( r ' ^)H( r ; t) represent the photon density in (r, i). These operators 
can be used to build the polarization correlation operators 

J ab = j d 3 rV;t(r,í)H(r,í) 

(18) 

dkÀl(k)À b (k), 
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where the last result follows immediately from Eq. (J17|) . The meaning of the matrix operator 
J = || J o6 || becomes clear when we calculate its trace: 

3 

= Jgg 



a=l 



(19) 



/ dk^2à\(k)à x (k) = N, 

J A=l 



where N denotes the photon-number operator. 

Now, by inserting the identity resolution Eq. (jïïj) in Eq. (JTHJ), we obtain 



J a b = J Jk^|0)(0|^4.fc + {multiparticle states} 



(20) 



dk|k,a)<k,6| (a, b =1,2, 3), 



v 



where the last equality holds only in the Hilbert spaces spanned by the one-photon states, 
V represents the set of detected modes, and we have defined the single-photon states |k, a) 



as: 



it|0) = |k,a), 



1,2,3). 



(21) 



From Eqs. (|19H2(J|) it is clear that Jn, J 2 2 and J 33 form a POVM (positive operator valued 



measure |4J) in the space spanned by the one-photon states. As we saw previously, the 
operators J a b's allow us to introduce a 3 x 3 correlation matrix operator 



/ j j f \ 

Jn J21 J31 



J 



\ 



(22) 



/ 



Jl2 J22 J32 
Jl3 J23 J33 

Now, for a light beam properly collimated around the direction = z, the 2x2 matrix 
obtained by extracting the first two rows and two columns from JJ = (J), coincides with the 
well known coherency matrix of the beam |5j, where the bracket average (•) is understood 
with respect to the state of the field. More generally, three independent 2x2 matrices can 
be extracted from J: 

a,b,c E {1,2,3}, 
c > ò, 
c, b 7^ a. 




(23) 
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In principle, each j( a ) = (j^) can be determined by measuring the Stokes parameters of 
the beam (either classical or quantum), with a polarization analyzer whose axis is parallel to 
e ( a ) p or example, for a beam propagating along the axis z, a set of generalized Hermitian 
Stokes operators jl|, |y] can be defined as 

S^Tria^J®}, (// = 0,...,3), (24) 

where the (p = 0, 1, 2, 3) are the normalized Pauli's matrices [7j. Then from Eqs. (j23H24j) 
it readily follows 

3 

J (3) = E^)' (25) 

where s M = (S^). Now we are ready to accomplish our initial task by introducing the 2x2 
effective reduced density matrix 

It is easy to check that when the set T> of the detected modes reduces to a single mode 
k || the defmition of p e ff above coincides with the well known polarization density 

matrix of a photon More generally, by using Eq. (jSJ, it is possible to introduce an 
effective 3 x 3_ reduced density matrix as p = J/Tr[J], which coincides with the one given by 
Peres et al. 



Single-photon scattering 

Let us consider now a genèric scattering process which transform the initial single-photon 
density operator p m in the output density operator p out . The most general linear transfor- 
mation between p m and p out , which leaves the density operator Hermitian and positive 
semidefmite, can be written 

p mt = J2paT a p íu T a \ (27) 

Aes 

where the scattering system has been represented by the ensemble S of scattering matrices 
{T A }, each of them occurring with probability pa > 0. If we insert Eq. (JjJJ into Eq. (|2Tj) 
we obtain 

1,2 r 

^ut = / dkdk'p^(k,k')J2pAr A \k,x)(k',\'\r A \ (28) 

A,A' J AeS 
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where 

1,2 



J2PAT A \k,X}(k' 1 X'\T A ^ = J2 í ^q^q , ^p A |q^)(q^|T A |k,A)(k' ) A'|T^| q ^)( q ^| 

Aes e,8' J Aes 



1,2 



E / dqdq'X;PA^(q,k)7;fí(k',q')|q6»>(q'(9'|. 



Aes 



(29) 

From the equation above, it is straightforward to see that we can write 

1,2 . 

p out = E / dqd q 'p$(q,q')|q0>(q , 0'l (30) 

where 

1,2 . 

Pet(<l, O = E^E / ^kdk%i(q, k)p* ,(k, k')T A ft (k', q'). (31) 
Ae<S A,A' ^ 

Since p*(k, k') = | \p\ x , (k, k') 1 1 (t = in,put) are 2x2 matrices, it is always possible to express 
them in the complete Pauli basis as p| 

3 

p t (k,k') = E^( k ' k ')^)> (32) 

At=0 

where we have introduced the two-mode Stokes parameters 5*(k, k') = Tr{o"( AÍ )p t (k, k')}. 
Since Eq. (j31j) can be written in matrix form as 



p out i 



Aes 

it is easy to write 

3 



q, q') = E^ / dkdk'T A (q, k)p in (k, k')T A t(k\ q'), (33) 



sTM) = É / rfkdk'^^^kOE^I^^^k)^)^^'.^)} 



i/=o " Aes 



3 

= £ / dkJk'Òk,k')E^m^(q,q';k,k') (34) 
u=o Aes 

3 

= E / rfkJk'M^(q,q';k,k / )5'Íf(k,k / ), 
where the four-mode density Mueller matrix 

m^(q, q'; k, k') = Tr K)T A (q, k)a H T A t( k ', q')} , (35) 
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is defined for a single ensemble realization A, while 



M^(q,q';k,k') = ^Pa <,(q, q'; k, k'), (36) 

represents the ensemble-averaged four-mode density Mueller matrix. It is easy to see that 
when the input state is a single-mode k state: 

,(k, k') = PA v5 (3) (k - k )5 (3) (k' - ko), (37) 

then 

5; n (k, k') = #*( 3 >(k - k )<^(k' - k ), (38) 

where p = \ \p\\>\ \ and s" 1 = Tr{o"( M )p}. In this case, the single-mode input and output Stokes 
parameters have the same functional relation as their classical counterparts: 

3 

S° ut (qo,qo) = / dkdk y M^(qD,q D ;k,k , )s^(k-ko)<j( 3 )(k , -ko) 

3 

= ^M^(q ,qo;k ,k )4 n (39) 



m 



U=0 
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